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ABSTRACT 

Inventory  systems  may  be  analyzed  using  techniques  which 
were  originally  devised  to  study  the  dynamic  response  of  servo- 
mechanisms.   This  paper  uses  these  methods  to  study  an  inventory 
system  using  exponential  smoothing  to  predict  a  demand  which  is 
assumed  to  have  a  linear  trend  with  uncorrelated  noise  components. 
The  statistical  properties  of  the  estimator  are  investigated,  and 
both  exact  and  limiting  expressions  for  its  mean  and  variance  are 
derived.   The  model  is  assumed  to  have  a  fixed  time  lag  between 
the  placing  of  a  replenishment  order  and  its  receipt,  and  the 
system  response  is  determined  for  impulse,  step  and  ramp  inputs. 
The  choice  of  a  smoothing  constant  is  discussed,  and  the  effect 
of  the  choice  on  the  variance  of  the  forecast  equation,  and  the 
response  time  of  the  system  is  analyzed. 
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1.  Introduction. 

This  paper  is  a  study  of  a  dynamic  inventory  system  using 
techniques  of  linear  systems  analysis(l).   The  model  is  an  inventory 
control  system  with  a  fixed  order  point,  and  a  known  fixed  time  lag 
between  the  placing  of  a  replenishment  order  and  its  receipt.   The 
demand  is  assumed  to  be  random,  with  a  linear  trend  and  uncorrelated 
noise  components.   Exponential  smoothing  is  used  to  forecast  future 
demand  (2).   The  reorder  rule  is  designed  to  minimize  inventory 
variation  for  any  sequence  of  forecasting  errors  (4).   The  controlled 
quantity  (inventory)  may  be  either  positive  or  negative.  A  negative 
inventory  represents  a  situation  where  the  demand  has  exceeded  the 
material  on  hand,  and  back  orders  have  been  placed  against  future 
deliveries.   In  order  to  facilitate  the  mathematical  manipulations, 
the  system  equations  will  be  transformed  using  generating  functions, 
and  the  response  of  the  system  to  deterministic  inputs  (impulse, 
step,  and  ramp)  will  be  investigated.   In  this  paper  the  term 
estimate  is  synonymous  with  predict. 

2.  Generating  Functions 

In  order  to  facilitate  analysis  of  the  system,  the  variables 
shall  be  treated  as  infinite  discrete  time  series,  and  their  gener- 
ating functions  will  be  obtained.  For  example  if  X  is  the  input 
(demand),  with  components  X(n) ;  n=l,2,...,  then  the  generating 
function  for  X  is  given  by, 
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An  important  characteristic  of  generating  functions  is  the 
following: 

If 

m 

(2-2)      Q-(/n)  -     Ji  fwhcm-&) 

where  g,  f  and  h  are  arbitrary  discrete  time  series  with  components 
g(n) ,  f  (n)  and  h(n)  then 

(2-3)      YL^^^     =    T   Z°     ffwhM) 
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or 


^T(2)    -  fT(i)  W(l) 


Thus,  the  convolution  operation  on  the  original  variables  is 
analagous  to  multiplication  of  the  transformed  variables. 

3.   Flow  Graphs 

Flow  graphs  are  compact  "visual  devices"  that  help  describe 
relationships  in  a  linear  system  where  the  basic  operation  performed 
is  convolution.   For  example  equation  (2-3)  may  be  represented 
graphically  as: 

i   o h    o  \ 


where  the  input  to  the  system  is  some  descrete  time  series  f,  and 
the  output  is  another  discrete  time  series  g  which  is  produced  by 
the  convolution  of  f  with  the  system  "transfer  function"  h.   The 
circles  at  the  input  and  the  output  of  the  system  are  called  "nodes" 
and  any  path  connecting  two  nodes  is  called  a  "branch".   The  purpose 
of  introducing  this  notation  is  that  more  complex  systems  may  be 
simplified  into  equivalent  systems  such  as  pictured  above  by 
applying  some  general  equivalence  rules.   For  example 


is  equivalent  to 


a+  4- 


because  of  linearity. 

A  network  which  is  similar  to  the  inventory  problem  is: 


where  i  is  the  input  to  the  network  and  o  is  the  output.   The 
letters  above  the  branches  represent  the  transfer  function  of  their 
respective  branches. 

This  is  similar  to  networks  found  in  electrical  control  systems 
and  methods  devised  for  analysis  of  these  physical  systems  work 
equally  well  in  the  analysis  of  the  inventory  system.  A  first  step 
in  the  analysis  is  to  find  any  "loops"  in  the  systems,  i.e.,  a  closed 
path  through  the  network  in  the  direction  of  the  arrows  which  touches 
no  node  more  than  once  in  one  traversal  of  the  "loop".   In  the  above 
example  cde  is  a  loop.   The  "loop  product"  is  defined  as  the  negative 
product  of  the  transfer  functions  of  the  branches  making  up  the 
loop.   In  the  above  example  -cde  is  the  "loop  product".   The  network 
determinant  (A)  is  defined  as  1  plus  the  sum  of  all  the  loop 
products  in  the  network.   In  this  case: 
A  =  1  -  c  d  e 

The  final  quantities  required  are  the  path  transmission  ?L   and 
the  path  determinant  A*.  .   P^  is  simply  the  product  of  the  branch 
transfer  functions  for  any  simple  path  between  the  input  node  and 
the  output  node  in  the  direction  of  the  branch  arrows.   In  the 
example  above  there  are  two  path  transmissions  P.  =  a  d 

P2  =  b  c  d 
The  path  determinant  is  1  plus  the  product  of  all  loops  in  the  net- 
work which  do  not  share  any  node  in  the  path.   In  the  above  example 

*1  "   A2  -  1 


The  network  transfer  function  Is  then  defined  as 


(3-D  T  =     ZL  P^  a^ 


A 

In  the  above  example 


(3-2)       X  =     q.c<  t  led 


This  result  will  be  used  in  the  analysis  of  the  inventory  problem. 

4.   Exponential  Smoothing 

The  operation  of  exponential  smoothing  is  a  weighted  average 
of  past  observations  used  to  predict  what  some  future  value  of  the 
variable  will  be.   The  weight  given  to  previous  observations 
decreases  geometrically  with  age.   The  forecast  equation  used  in 
this  paper  is  a  linear  combination  of  single  and  double  exponential 
smoothing.   It  will  be  shown  that  a  linear  combination  of  single 
and  double  smoothing  yields  an  asymtotically  unbiased  estimate  of 
the  demand  when  the  demand  is  assumed  to  have  a  linear  trend. 

The  Single  Exponential  Smoothing  operator  is  given 
by: 


=  <X(t)   +-  3  L<X(t-o  +  6  ^  (*>J 
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Where  X0  Is  the  initial  estimate  of  the  demand  and   p  s   |  —  «c 
ind  has  no  trend  component,  i.e.; 


yet)  =  a,  1 6         t  >o 


Where  a  Is  a  constant  and  c  is  random  with 


E  Ut]  -o       y  t 
EUt'J=<$'    vt 


ElMil  s  o    V-issj 


Then, 


t-» 


=    a  -    a/s*    +   (3*e  Lx.] 

If  the  initial  estimate  of  the  demand  is  an  unbiased  estimate 

based  upon  past  observations  then 

E  [S(K)]  s   3.    V  tl  and,  in  any  event, 

JU.  E  L  V^  =  «- 


Since  smoothing  is  a  linear  operation: 

(4-3)   Vx;  =  |c*u   +  St(«) 

*-' 


<=o 


r    a-  +  Stu) 


(4-4) 


Sa  (x)  =  a2  +  2a.  Sen  +  52(6) 

t  *        t 

1  t_,     K-e>      t-K       ^=,  Js0    t-^t-j 
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Where  O^  is  the  variance  of  the  initial  estimate  of  demand. 


Hence, 


V*tlstc«]  ■     F[St2(x)]  -  [E(st(«)|l 


It  follows  that: 


The  above  results  indicate  that  exponential  smoothing  gives  an 
asymtotically  unbiased  estimate  of  the  demand  and,  if  the  criterion 
for  choosing  the  smoothing  constant  were  minimum  variance,  the  choice 
«*  =  0  would  be  optimal.   This  of  course  would  be  the  same  as  using  the 
initial  estimate  of  demand  for  all  future  estimates  and  would  negate 
the  most  desirable  trait  of  exponential  smoothing,  i.e.:  adjusting 
the  estimate  when  the  time  demand  changes.   It  does  however  indicate 
that  where  there  is  a  high  degree  of  confidence  in  the  initial  estimate, 
and  it  appears  that  changes  in  demand  are  unlikely,  a  small  value 
of  the  smoothing  constant  is  preferable. 

5.   Double  Exponential  Smoothing 

For  inventory  situations  where  the  demand  is  assumed  to  have 
the  form: 

X(t)  "Q-  +  c.  as  discussed  above,  the  single  smoothing 
operator  gives  a  suitable  estimate  of  future  demand.   If  demand  is 
assumed  to  have  a  linear  trend  component,  then  it  will  be  shown 
that  both  single  and  double  smoothing  give  biased  estimates  of  the 
demand  even  in  the  limit,  but  a  linear  combination  of  the  two 
operators  yields  an  asymtotically  unbiased  estimate. 
For  the  linear  model: 

X(t)  s   a  ♦  bt +£  where  a  and  b  are  constants  and, 

Etcj  =  o 

EL^]  --  <rft*  Vt     with 


EtMil*  °    V>L*;j. 


Let  f  :  a  +  bt  so  that, 

St(X)  *  St  (f)  ^St  (ft)  by  linearity. 

Now,     ft.K=      <K  +  lrt~  JrK     -    f^-** 

si««    E[S(o]«    El^X'eV      +fi\]  •  o 

(5-2)    £  L  V^J  •   £  [  V  ?')  *  £  L^>]  --  V  j) 

(5=3)       £*t[E[V*)]-    £^J     *    X^fV?)-^-^] 

=    -JUL 

The  above  result  shows  that  the  single  smoothing  operator  does 
not  give  an  unbiased  estimate  of  the  demand  in  the  linear  model 
even  asymtotically.  Moreover,  after  a  sufficient  time  period  to 
allow  the  exponential  terms  to  become  negligible  the  expected  value 
of  the  single  smoothing  estimator  may  be  written  as; 

<5-4>      ECStxjJ  =    CL+it-ii. 


Similarly  for  double  smoothing: 


(5-5)    5(X)  =   5  ({)  f  J+  <o 

t         t   •       t 


0 
K=   °  Ki  *  fe  S  ft 


? 


2         f  r  /6 


Sine.   EtS^«]  =  0 

(5-7)    E[StW(x)J-  5tl*(T) 

ana  ^E[^M^  ^tflH-tJ-^ 

Here  again  the  estimate  is  biased  and  the  limiting  mean  may  be 
written  approximately  as: 


(5 


8)     tLS^cx;]   =  O^  +  Jrt-  ibg. 
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Using  these  approximations  it  can  be  seen  that: 


(5-9)  a.    '-     £    J)  [xi     -     S      Cxi 

(5-io)  >     =   ~    I  S  tx]  -   S  C2cxi  1 

A     L    t  t        J 


are  asymtotically  unbiased  estimates  of  the  parameters  in  the  linear  -T 
model.   Although  the  estimates  are  not  unbiased  the  bias  does 
decrease  exponentially  with  time.   Using  these  estimators,  the  fore- 
cast of  demand  at  time  (t  +V)  becomes: 

(5-u)      stc x j r)  .    [Z  +  y]  stoo  -  L -  +  ^J  s*'c*> 

The  generating  functions  for  exponential  smoothing  will  now 
be  introduced.  As  given  in  equation  (4-1),  single  smoothing  may 
be  written  as: 

t-l 


K*o 


which  is  the  convolution  of  the  single  smoothing  operator  with  the 
discrete  time  series  of  demand.   In  the  limiting  case,  i.e., 
where  the  initial  demand  time  becomes  negligible  (or  setting  X  x  0) , 
the  transformed  equation  may  be  written  as: 


Su>  =  Xrc*)  STc«) 


ii 


where  S   (Z)  is  the  transform  of  an  exponential  decay  function,  i.e. 


may  be  transformed  to 


00 

(5-12)     5(a)  -  <  l_   £  i       =   ,  . 

Similarly  for  second  order  smoothing: 

may  be  recognized  as  the  convolution  of  first  order  smoothing  with 
itself  for  large  t,  and  in  the  transform  domain  this  is  the  square 
of  the  transform  for  first  order  smoothing: 

It  now  remains  to  find  the  variance  of  the  forecast  equation. 
First  the  variance  of  the  second  order  smoothing  operator  will  be 
found. 

For  the  linear  model: 

5.  ex)  --  S  m  4    v\i) 

Hence  Var  L^CX)}  *  £  [SW(fc)f   sinc«  S^2)  (  )  is 
and   E  t  S^(€)]  ■  O 


constant 
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Now, 

where  £. 


Sl\k) --   oc  zi  ys Ks  co    1-3  *J 


•  r  til 


•      =        °o    ^ 


and,  S    CO   .    *£     £*«  +     £ 

So   that,        ^ 

t 


K'°  ^  *"*"*  fro         •        ' 


The  second  moment  Is: 

(5. '  '*-     '***■ 


where  tf~  is  the  variance  of  the  initial  estimate  for  first  order 
smoothing  and  <r.  is  the  variance  for  second  order  smoothing. 
The  summation  in  the  first  term  of  equation  (5-14)  may  be  evaluated 
directly  or  it  can  be  recognized  that  this  is  the  sum  of  the  first 
t  squared  terms  of  the  impulse  response  of  the  second  order  smoothing 
operator,  i.e.; 

cT»>      -  •<•* 
5  (,)  -  Ci-G*)x 

where  S  \t)   is  the  Z  transform  for  the  second  order  smoothing 
operator.   If  we  denote  the  impulse  response  by  F(z); 


r<«  =  x7«  $% 
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where  X  (Z)  r  1  (the  Z  transform  of  a  unit  impulse  at  t  ■  0) . 


The  impulse  response  is: 


and  the  sum  of  the  first  t  squared  terms  is 
*"'  u    *"'      *   2  it 


vhere  T  =/?>  ;  therefore:  *r° 

\l         rt2t   *      2^2* 


(5-15) 


n^«r^-;+^vv  +  -v[T^ 


*<'■">       „j,       Ut-ti)  _zt  i 


iW 

and  in  the  limit  this  becomes: 


L'-^3  U-fJ1] 


u.    ftSlW=     •<  C '  +  P*l      a-* 

so   that 

(5-16)   Var    [  *  WU>]    "-       L    L  S^Ci)]  ' 

The  transfer  function  approach  may  also  be  used  to  evaluate  the 
variance  of  the  forecast  equation.   In  this  case 

L         (i  J  '-(Si        L        $  J(,-(s^)i 
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The  impulse  response  is: 


and  the  sum  of  the  squares  of  the  response  terms  is 

The  limiting  or  asymtotic  variance  of  the  estimator  then  becomes 

(5-17)  var  f  S  M>|  =        l*rt  +  ™x+n«*Wr+iSrx}t9:* 

It  is  then  apparent  that  the  proper  linear  combination  of  single 
and  double  exponential  smoothing  will  give  an  asymtotically  un- 
biased estimate  of  demand  where  there  is  a  linear  trend,  and  the 
noise  components  are  uncorrelated.  As  can  be  seen  in  the  derivation 
of  the  moments  of  the  smoothing  operator  the  variance  of  the  operator 
is  a  weighted  sum  of  the  variance  of  the  observed  demand  and  the 
variance  of  the  initial  estimate.  As  more  observations  are  obtained 
the  variance  then  approaches  the  limiting  variance  as  given  by 
equation  (5-17)  for  the  forecast  equation  of  the  linear  model.   In 
the  limiting  case  it  can  be  seen  that  small  values  of  U.   give  the 
smallest  variance.   It  will  be  seen  later  in  the  paper  that 
although  small  variance  is  a  desirable  characteristic  of  any  fore- 
cast, the  response  of  the  forecast  equation  to  changes  in  demand 
pattern  is  also  a  function  of  u.   and  if  u.   is  chosen  to  be  small  the 
model  will  be  slow  in  responding  to  changes  in  the  demand. 
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6.  Fixed  Lag  Inventory  Model: 

The  problem  to  be  considered  is  an  inventory  system  where  there 
are  fixed  order  points.   In  practice  this  would  correspond  to  an 
activity  which  finds  it  more  economical  or  necessary  to  adhere 
to  a  regular  order  schedule,  and  can  only  place  replenishment 
orders  at  the  beginning  of  each  input  period.   Specific  character- 
istics of  the  system  are: 

1)  The  item  has  high  usage,  and  there  are  indications  of  a  linear 
trend  in  the  demand. 

2)  There  is  no  seasonal  variation  in  demand,  and  no  correlation 
in  the  error  terms. 

3)  The  order  points  are  fixed  and  equally  spaced. 

4)  When  demand  exceeds  the  inventory  on  hand,  back  orders  are 
placed  against  future  deliveries. 

5)  There  is  a  fixed  time  lag  between  placing  a  replenishment  order 
and  receipt  of  the  order.  Time  lag  is  measured  in  demand  periods. 

6)  Time  is  measured  from  the  end  of  the  demand  period. 

7)  A  replenishment  order  placed  at  time  t,  may  be  used  to  fill  demand 
generated  during  demand  period  (t,  +  T  +  1)  ,  where  the  time  lag  is  T. 

8)  Oversupply  may  be  returned  to  the  vendor;  i.e.,  a  negative  order 
quantity  may  be  generated. 

9)  The  forecast  equation  introduced  in  section  (5)  will  be  used 
to  predict  future  demand. 

10)  The  forecast  at  time  k  will  predict  demand  during  period  k  +  T  *  1, 
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The  following  equations  describe  the  model: 


(6-1) 

(6-2) 

where: 
(6-3) 


i(t)  =  i.  +  Lrl  -  ixci) 


<--  I 


5t(x;r;=  l**£]5cx)  -  [,  +  <f]  s^cxi 


S  co    =    «<.  X(t)    +  (3  S  (x) 


(6-4) 


(6-5) 


(6-6) 


S    (X)    ~      •<•  Sex)    +    6  S  2cx) 

Oct)   =      2T  s  ex;       -  foit-d)-  let) 
RCO=     O(t-CT-n)) 


The  flow  graph  for  this  model   is: 
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where: 
(6-7) 


^  [(2+^T)-  (Z+*(T-|))£l 


which  is  a  linear  combination  of  the  transforms  for  single  and 
double  exponential  smoothing.   The  network  determinant  is: 


A--     i  + 


I 


■-M 


The  path  transmissions  and  path  determinants  are: 

i  -  %  x 

and  the  system  transfer  function  is 

(6-8)  S--   P.A.  +  P^x  .  ^T^-*TT[(^T;-fc+^T-,>)*]    ,_£*' 

The  purpose  of  describing  the  inventory  model  in  terms  of 
flow  graphs  and  generating  functions  is  to  determine  the  dynamic 
response  of  the  system  to  possible  demand  patterns,.   It  was  seen 
that  in  examining  the  moments  of  the  forecast  equation  the  variance 
of  the  forecast  equation  was  a  function  of  the  smoothing  constant 
*C  and  the  lead  time  T.   On  examining  the  dynamic  response  of  the 
system  to  deterministic  inputs  it  will  be  seen  that  the  same  para- 

J 
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meters  determine  the  response  time  (time  until  all  back  orders  are 
filled).   The  model  will  be  tested  with  three  deterministic  inputs 
(impulse,  step,  ramp)  and  the  model  will  be  checked  for  response  to 
these  inputs,  as  well  as  the  steady  state  error. 

In  order  to  facilitate  the  determination  of  steady  state  error 
the  following  theorem  from  (3)  is  introduced. 

Let  F  Be  the  transform  of  f  and  suppose 

F(Z)  converges  as  Z  -»  1  where  Z  is  real  and 

1Z1<  1  and  f(n)  =  0  for  n<0.   Then: 

lim  f(n)  =  f(*>)  =  lim  (1  -  Z)  F(z) 
h-»«»  2-»| 


6.1   Impulse  Response 

The  impulse  response  for  the  inventory  model  in  transform 
notation  becomes 

Hence  by  the  above  theorem; 

I(05)=  JL^,  /I-*)   \(/-2*)*  Z^'lil+dT  )~    U  +  *[T-H)  g 


Here  I  (<»)  denotes  the  "steady  state"  inventory. 

In  order  to  obtain  the  inverse  transform,  the  transfer 
function  may  be  expanded  into  a  sum  of  partial  fractions: 

-£_   .  *[n-i-r]    _L__  1    /-  fc7"*1 


TT,        2T+'    aT4l  f  -P      .   otLl  +  ^1    .  _J ] 


i  -  % 
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Then  solving  for  the   inverse: 

ICt)    =      -/3LP  H   (t-CT+O)   ~    /3  MCt-iCT+.)J 


f   ^[.-ht]   j^  ct  -r)fZ(t-7)-(t-  Cir+Oj^Iut^iTfo] 


(6-9)  *   ^  [1+JTi   r  . .       ^ti-r 

where     ju.(K)  -  |       K  >  Q 
-  O       K  <o 
In  addition  to  the  fact  that  there  is  no  steady  state  error, 
the  expression  for  the  impulse  response  also  shows  that  if  the  system 
was  receiving  a  constant  demand  for  a  long  period  of  time  and  there 
was  a  one  time  increase  in  that  demand  the  system  would  respond 
in  the  minimum  amount  of  time  since: 

Im  »   -  i     ;  t  4.  t+i 

Kr+o  =    *c  (z  +  jlt) 

In  other  words,  there  would  be  a  back  order  for  a  period  equal 
to  the  delivery  time  but  all  demand  would  be  met  by  time  T  +  1, 
and  then  in  fact  there  would  be  an  excess  of  material  (  *c  £  2  +  <*  T/y 


6.2  Response  T®  A  Step  Imiput  (  X  ( £)   -  yz~£  ) 

Testing  the  system  with  a  step  input  gives  an  indication 
of  system  performance  when  there  is  a  sudden  increase  in  demand 
and  demand  continues  at  this  new  higher  level.   The  results  which 
are  derived  in  this  example  apply  equally  to  the  case  where  there 
is  a  sudden  drop  in  the  demand  level. 
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I  f«  -   S  trm   -    4 


Using  the  same  partial  fraction  expansion  as  above  the  inverse 
transform  yields: 


(6-10) 


t-C7~*-i)  t-z(T+/) 

lit)  '-   --£[(!-£)  M-Ct-iTto)  -  (t-  P    )  M.Ct-zCn-1)) 


U. 


+  Liiilii/-  [i+  (t-TV]|g  ~^u-c7>,;;-  [,  - 

r  t-27-/  n 

[  '  +  Ct-H-  z(t>/)*]0   AtCt-  Z<T-m))J 

+  zCt-T)^(t-cr^o)-  (t-  (2_T+/))xCt-  2(r^/)) 

As  expected  there  is  again  negative  inventory  after  the 
initial  change  in  demand.  At  the  start: 

Tct)  =  -  (t-n)  ;  Oit  <Tfi 

The  first  order  is  received  at  (T  +  1)  in  the  amount  2  +  «t(l+<T). 
It  appears  then  that  if  the  criterion  for  selecting  a  smoothing 
constant  was  minimum  response  time,  <<  s  1  would  be  the  optimal 
choice.  This  would  be  the  extreme  case  where  no  weight  was  given 
to  past  demand  and  the  forecast  was  based  upon  only  the  last  recorded 
demand.  This  of  course  would  be  impractical  since  actual  demand 
will  normally  have  a  noise  component  (i.e.  be  probabalistic  rather 
than  deterministic)  and  as  shown  previously  a  choice  of  ii   -    1 
would  make  the  variance  of  the  forecast  equation  maximum. 
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It  Is  then  evident  that  no  single  value  of  the  smoothing  constant 
will  satisfy  both  the  requirement  of  minimum  forecast  variance  and 
minimum  response  time. 

6.3  Response  T©  A  EUmp 

This  is  the  type  input  for  which  the  forecast  equation 
was  designed.  The  reasons  for  looking  at  the  step  and  Impulse 
responses  were  to  see  how  the  system  performed  under  other  than 
designed  conditions.  Since  the  system  is  linear  and  time  in- 
variant the  performance  of  the  system  to  any  combination  of 
inputs  may  be  analyzed  by  looking  at  each  input  separately  and 
obtaining  the  actual  response  by  adding  the  system  response  to  each 
of  the  components  of  the  input. 

For  the  ramp  input:   X  ( £)  —  i    — \t- 
Expanding  X  (*)  as  before  and  solving  for  J  ft; 


**T+i  K=o 


t-2T-l 
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It  is  difficult  to  analyze  response  time  from  the  above 
expression  but  if  a  ramp  input  is  looked  upon  as  a  sequence  of 
impulse  inputs  where  the  magnitude  of  the  impulse  at  time  t  is 
(t  -r  1)  the  superposition  property  of  a  linear  time  invariant 
system  leads  to  the  conclusion  that  the  larger  the  value  of  oC  , 
the  smaller  the  response  time. 

6.4  Order  Rule 

As  shown  by  Vassian  (3),  the  order  rule  utilized  in  the 
present  model  gives  minimum  inventory  variation  for  any  sequence 
of  forecasting  errors. 

7.   Distributed  Lag  Inventory  Model 

The  preceding  model  with  a  fixed  lag  time  was  characterized 
by  a  system  response  which  was  a  function  of  the  smoothing  constant 
and  an  inventory  which  tended  to  overshoot  the  demand  when  <<   was 
chosen  to  be  large  in  order  to  decrease  the  response  time.   The 
same  problem  will  now  be  examined  when  the  lag  time  is  distributed 
rather  than  fixed  (5).   In  the  context  of  the  inventory  model  we 
now  have  that; 

(7-D         R<t)  =•    £ot*)  frt-i) 

where;       f  (t  j  \>)   *  p  <^ 

is  the  probability  that  an  order  is  received  exactly  t  days  after 
it  is  placed. 
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The  order  rule  will  now  be; 


(7 


-2)      Oct)  -    5  (*;t>  -  let; 


and  the  forecast  equation  will  predict  the  demand  at  T  «   '  I  V 
(the  expected  delivery  time). 
The  flow  graph  Is; 


XT<*) 


-*Q 


15(*) 


-K>- 


fk> 


i-£ 


1T(*) 


and 


r< 


*)    - 
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This  network  has  the  transfer  function: 


where  r.  and  r«  are  the  roots  of  the  second  degree  polynominal  in 
the  denominator  of  S,  and; 


a, 


\ 


fHf^-^J2  IV,] 

[Hf^-/^[-^] 


,  where     IT  =  /  5" V  -  ^  9  +-  / 
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7.1     Impulse  Response 

IT(*>  =  S 


(7-4) 


let)-   -±l>\     &«*,)?-  A.^*,*™ 


P: 


I  <<*>)  =  o 

Just  as  with  the  fixed  time  lag  model  there  Is  no  steady  state 
error,  but  the  transient  response  does  exhibit  a  property  not 
evident  in  the  fixed  lag  model.   For  0  £  <b  C  .2.  R  is  real  and  all 
of  the  terms  in  J  (t)  are  exponential  decay.   For  .2^0.  <  \,0      R  is 
imaginary  and  the  coefficients  A«  and  A,  are  related  byj 


/v» 


A3  =  Af 

and         Y",  "  ^t 

so  that  when  written  in  exponential  form: 

A,»    Ce**  Af-     Ce'}& 

and  the  last  two  terms  in  I{t)  became 

J(t)- CD  [c*        *  c     J 


-CDclL^+  Ct+04>] 
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The  response  then  is  the  sum  of  terms  containing  exponential 
decay  and  damped  oscillation. 

The  results  for  a  step  and  ramp  input  as  seen  below,  also 
contain  terms  which  are  damped  oscillations  for  the  values  of  q 
given  above.   This  result  is  not  unexpected.   Damped  oscillation 
is  present  in  physical  systems  in  which  there  is  delay  in  response. 
The  lag  in  receipt  in  the  inventory  system  is  analagous  to  inertia 
in  a  mechanical  system  or  inductance  in  an  electrical  system,  and 
both  of  these  physical  systems  will  have  transient  responses 
containing  damped  oscillation  under  certain  conditions  which  are 
analagous  to  the  relationship  of  q  in  the  inventory  system. 

The  steady  state  error  for  both  the  step  and  ramp  input  is 
zero  in  the  distributed  lag  system  and,  other  than  this  tendency 
to  oscillate,  the  results  are  similar  to  those  in  the  fixed  lag  model. 

7.2  Inventory  Function  For  A  Step  Input  (Distributed  lag); 
when  demand  is  a  step  the  response  is; 

(7-5)      i(t)  -  jM/W-fiv  [ i-o +  «♦.>.<] s**'i 

7.3  Inventory  Function  For  A  Ramp  Input  (Distributed  lag); 
when  demand  is  a  ramp  the  response  is; 
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;,6)   I(t)s   -A  [t-£[H»*]j*j*[t-  *♦»♦*]*] 

a  f  r,      fr  +  i  _   '  .      r 


-t 


7.4  Summary 

The  distributed  lag  model  is  probably  a  more  realistic 
description  of  a  practical  inventory  problem  in  that  the  delivery 
time  for  replenishment  orders  would  seldom  be  a  fixed  constant 
as  was  assumed  in  the  first  model.   Practical  inventory  problems 
are  characterized  by  the  damped  oscillations  of  the  distributed 
lag  model.  When  the  delivery  time  is  not  fixed,  it  should  be 
expected  that  there  will  be  periods  of  shortage  and  inventory 
excess  caused  by  variations  in  delivery  time.  As  can  be  seen  in 
the  above  equations  response  time  is  now  a  function  of  the  smoothing 
constant^  and  |»  s  the  parameter  of  th<?  geometric   distribution  of 
delivery  time.   The  inventory  equations  could  be  analyzed  directly 
in  order  to  choose  an  acceptable  value  for  the  smoothing  constant, 
but  as  long  as  it  is  recognized  that  there  will  be  some  oscillation 
no  matter  how  the  smoothing  constant  is  chosen  it  does  not  appear 
to  be  invalid  to  use  the  fixed  lag  formulation  of  the  problem,  where 
the  fixed  lag  is  chosen  to  be  the  expected  value  of  the  random 
variable  (lag  time). 
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8.   Selection  of  a  Smoothing  Constant 

In  choosing  a  smoothing  constant  for  the  model  it  must  be 
recognized  that  no  single  value  of*  is  optimal  for  all  demand 
patterns.   The  impulse  response  for  the  fixed  lag  model  shows  the 
relationship  which  exists  between*  and  the  time  lag.   It  will  be 
recalled  that  for  this  model: 

l(t)   *  -|    0£t   <   T+l 

X(T+o  =     <(*  *  *t; 

This  indicates  that  when  there  is  noise  present  in  the  demand  a 
long  time  lag  will  tend  to  amplify  the  fluctuations  in  inventory 
caused  by  the  noise.   This  can  be  offset  somewhat  by  choosing  a 
smaller  «<  when  T  is  increased,  but  when  «<  is  reduced  the  response 
time  is  increased.   The  effects  of  decreasing  ^  however  would  be 
felt  in  the  sensitivity  of  the  forecast  to  changes  in  parameter. 
If  the  slope  of  the  trend  component  were  to  increase,  a  smaller 
value  of  <   would  cause  the  forecast  to  lag  behind  the  demand  for 
a  longer  period  of  time.  Any  final  choice  of  •<  would  also  have 
to  take  into  consideration  the  "cost"  of  stock  outages  vs  the 
holding  cost  of  excess  inventory.   The  model  considered  controlled 
the  inventory  about  the  zero  level  and  in  most  practical  cases  a 
non  zero  safety  level  of  stock  would  be  maintained.  Any  non  zero 
safety  level  would  alleviate  the  response  time  problem  since  this 
buffer  stock  would  compensate  for  forecast  errors  during  a  period 
of  change  in  demand  pattern.  All  of  these  factors  seem  to  favor 
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a  small  smoothing  constant  to  increase  system  stability,  and  a 
safety  stock  level  to  compensate  for  the  increased  response  time. 

9.   Conclusions 

The  use  of  flow  graphs  and  generating  functions  provide  a 
simple  method  of  obtaining  some  insight  into  the  characteristics 
of  an  actual  inventory  system.   Testing  the  model  with  deterministic 
inputs  allows  the  inventory  manager  to  check  his  reorder  rule  for 
stability  and  to  decide  if  his  forecast  equation  is  satisfactory 
as  regards  response  to  fluctuations  in  demand  patterns.   The  primary 
advantage  of  using  exponential  smoothing  to  predict  demand  is  its 
flexibility  and  its  economy  of  computer  storage  space  where  it  is 
used  in  a  mechanized  activity.   The  flexibility  allows  the  inventory 
manager  to  make  trade  offs  between  response  time  and  forecast  var- 
iance by  simply  adjusting  the  smoothing  constant.   If  the  demand 
pattern  changes  at  any  time,  he  may  increase  or  decrease  oL   ,  and 
if  this  a  temporary  change  in  pattern  he  may  at  any  time  return  to 
the  former  value  of  •<•  .   The  economy  of  storage  space  is  apparent 
when  exponential  smoothing  is  compared  to  the  use  of  a  moving 
average  to  predict  demand.  Where  the  moving  average  estimator 
using  the  last  N  observations  of  demand  would  require  N  storage 
locations  for  each  item  of  inventory,  exponential  smoothing  requires 
only  one  location  to  retain  the  last  smoothed  value.   In  this 
regard,  if  a  moving  average  estimator  provides  satisfactory  pre- 
dictions for  a  given  inventory  situation,  an  exponential  smoothing 
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estimator  can  easily  be  constructed  which  will  have  the  same 
limiting  variance  as  the  moving  average  estimator. 

For  example  if  a  moving  average  has  been  used  to  predict  demand 
for  a  constant  model;  i.e.: 

where  the  noise  is  uncorrelated  with  variance  <^  .   The  moving  average 
estimator  is 


and 


(8-2)    Y**  [  7  )  -  JSl 


as  shown  in  section  (4)  the  limiting  variance  for  single  smoothing 
in  the  case  of  the  constant  model  was: 


var  lSt(«J  =   77^"°** 


so  that  if  *  is  chosen  so  that 

(8-3)     ■«*  "  Z/N+l 
the  limiting  variances  of  the  two  estimators  are  identical.  A 
similar  criterion  could  be  established  for  the  linear  model  and  «*. 
could  be  chosen  to  make  the  variance  of  the  estimate  approximately 
the  same  as  that  for  an  acceptable  moving  average  operator. 
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